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Abstract — This article provides an overview of tlie state-of- 
art results on communication resource allocation over space, 
time, and frequency for emerging cognitive radio (CR) wire- 
less networks. Focusing on the interference-power/interference- 
temperature (IT) constraint approach for CRs to protect pri- 
mary radio transmissions, many new and challenging problems 
regarding the design of CR systems are formulated, and some 
of the corresponding solutions are shown to be obtainable by 
restructuring some classic results known for traditional (non-CR) 
wireless networks. It is demonstrated that convex optimization 
plays an essential role in solving these problems, in a both 
rigorous and efficient way. Promising research directions on 
interference management for CR and other related multiuser 
communication systems are discussed. 



I. Introduction 

In recent years, cognitive radio (CR) networks, where CRs 
or the so-called secondary users (SUs) communicate over 
certain bandwidth originally allocated to a primary network, 
have drawn great research interests in the academic, industrial, 
and regulation communities. Accordingly, there is now a 
rapidly growing awareness that CR technology will play an 
essential role in enabling dynamic spectrum access for the next 
generation wireless communications, which could hopefully 
resolve the spectrum scarcity vs. under-utilization dilemma 
caused by the current static spectrum management polices. 
Specifically, the users in the primary network, or the so- 
called primary users (PUs), could be licensed users, who have 
the absolute right to access their spectrum bands, and yet 
would be willing to share the spectrum with the unlicensed 
SUs. Alternatively, both the PUs and SUs could equally 
coexist in an unlicensed band, where the PUs are regarded 
as existing active communication links while the SUs are new 
links to be added. A unique feature of CRs is that they are 
able to identify and acquire useful environmental information 
(cognition) across the primary and secondary networks, and 
thereby adapt their transmit strategies to achieve the best 
performance while maintaining a required quality of service 
(QoS) for each coexisting active primary link. Depending 
on the type of cognitive knowledge collected (e.g., on/off 
statuses of primary links, PU messages, interference power 
levels at PU receivers, or primary link performance margins) 
and the primary/secondary network models of interests (e.g., 
infrastructure-based vs. ad hoc), many new and challenging 
problems on the design of CR networks can be formulated, as 
will be reviewed in this article. 

To date, quite a few operation models have been proposed 
for CRs; however, there is no consensus yet on the terminology 
used for the associated definitions H], 121, Generally 
speaking, there are two basic operation models for CRs: 



Opportunistic Spectrum Access (OSA) vs. Spectrum Sharing 
(SS). In the OSA model, the SUs are allowed to transmit over 
the band of interest when all the PUs are not transmitting 
at this band. One essential enabling technique for OSA- 
based CRs is spectrum sensing, where the CRs individually 
or collaboratively detect active PU transmissions over the 
band, and decide to transmit if the sensing results indicate 
that all the PU transmitters are inactive at this band with 
a high probability. Spectrum sensing is now a very active 
area for research; the interested readers may refer to, e.g.. 
El, 13, l6l. El for an overview of the state-of-art results in 
this area. As a counterpart, the SS model allows the SUs to 
transmit simultaneously with PUs at the same band even if 
they are active, provided that the SUs know how to control 
their resultant interference at the PU receivers such that the 
performance degradation of each active primary Unk is within 
a tolerable margin. Thus, OSA and SS can be regarded as the 
primary-transmitter-centric and primary-receiver-centric dy- 
namic spectrum access techniques, respectively. Consequently, 
there will be an inevitable debate on which operation model, 
OSA or SS, is better to deploy CRs in practical systems; 
however, a rigorous comparative study for these two models, in 
terms of spectrum efficiency and implementation complexity 
tradeoffs, is still open. Generally speaking, SS utilizes the 
spectrum more efficiently than OSA, since the former supports 
concurrent PU and SU transmissions over the same band 
while the latter only allows orthogonal transmissions between 
them. Moreover, the receiver-centric approach for SS is more 
effective for CRs to manage the interference to the PU links 
than the transmitter-centric approach for OSA. 

Hence, the SS model for CRs will be focused in this article. 
It is worth noting that the optimal design approach for SS- 
based CR networks should treat all coexisting PU and SU 
links as a giant interference network and jointly optimize their 
transmissions to maximize the SU network throughput with a 
prescribed PU network throughput guarantee. From this view- 
point, recent advances in network information theory [8 1 have 
provided promising guidelines to approach the fundamental 
limits of such networks. However, from a practical viewpoint, 
the centralized design approach for PU and SU networks is not 
desirable, since PU and SU systems usually belong to different 
operators and thus it is difficult, if not infeasible, for them to 
cooperate. Consequently, a decentralized design approach is 
more favorable, where the PU network is designed without the 
awareness of the existence of the SU network, while the SU 
network is designed with only partial knowledge (cognition) 
of the PU network. 

Following this (simplified) decentralized approach, there are 
furthermore two design paradigms proposed for SS -based CRs. 
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One is based on the "cognitive relay" concept Q, where 
the SU transmitter allocates only part of its power to deliver 
the SU messages, and uses the remaining power to relay 
the PU messages so as to compensate for the additional SU 
interference experienced at the PU receiver However, this 
technique requires non-causal knowledge of the PU messages 
at the SU transmitter, which may be difficult to realize in 
practice. In contrast, a more feasible SS design for the SU 
to protect the PU is to impose a constraint on the maximum 
SU interference power at the PU receiver, also known as the 
"interference temperature (IT)" constraint ifTOl . by assuming 
that the SU-to-PU channels are either perfectly known at the 
SU transmitters, or can be practically estimated. 

In this article, we will focus our study on the IT-based SS 
model for CRs, namely the IT-SS, as it is a more feasible 
approach compared with other existing ones. In a wireless 
communication environment, channels are usually subject to 
space-time-frequency variation due to multipath propagation, 
mobility, and location-dependent shadowing. Thus, dynamic 
resource allocation (DRA) becomes an essential technique for 
CRs to optimally deploy their transmit strategies to maximize 
the secondary network throughput, where the transmit power, 
bit-rate, bandwidth, and antenna beam should be dynamically 
allocated based upon the available channel state information 
(CSI) of the primary and secondary networks. In particular, 
this article will focus on DRA problem formulations unique 
to CR systems under the IT-SS model, and the associated solu- 
tions that are non-obvious in comparison with existing results 
|[m known for the traditional (non-CR) wireless networks. 
More importantly, we will emphasize the key role of various 
convex optimization techniques in solving these problems. 

The remainder of this article is organized as follows. Section 
inipresents different models of the SU network coexisting with 
the PU network, and various forms of transmit power and 
interference power (IT) constraints over the SU transmissions. 
Section |III| is devoted to the spatial-domain transmit optimiza- 
tion at the SUs for different SU networks subject to transmit 
and interference power constraints. Section |IV] extends the 
results to the more general case of joint space-time-frequency 
transmit optimization of the SUs, and addresses the important 
issue on how to optimally set the IT thresholds in CR systems 
to achieve the best spectrum sharing performance. Finally, 
conclusions are drawn and future research directions are 
discussed in Section IV] 

Notation: Lower-case and upper-case bold letters denote 
vectors and matrices, respectively. Rank(-), Tr(-), | ■ |, (•)^^, 
(•)^, and (•)^/^ denote the rank of a matrix, trace, determinant, 
inverse, Hermitian transpose, and square-root, respectively. 
I and denote an identity matrix and an all-zero matrix, 
respectively. Diag(a) denotes a diagonal matrix with diagonal 
elements given in a. E(-) denotes the statistical expectation. 
A circularly symmetric complex Gaussian (CSCG) distributed 
random vector with zero mean and covariance matrix S is 
denoted by CJ\f{0,S). C™^" denotes the space of m x 71 
complex matrices. || ■ || denotes the 2-norm of a complex 
vector. Rc(-) and Im(-) denote the real and imaginary parts 
of a complex number, respectively. The base of the logarithm 
function log(-) is 2 by default. 
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Fig. 1. CR networks: (a) infrastructure-based; (b) ad hoc. 

II. CR Network Models 

We consider two general types of CR networks, which are of 
both theoretical and practical interests: One is infrastructure- 
based, as shown in Fig. ^a), where multiple secondary ter- 
minals communicate with a common secondary node referred 
to as the secondary base station (S-BS); the other is ad hoc, 
as shown in Fig. [TJb), which consists of multiple distributed 
secondary links. In both types of CR networks, there are 
coexisting primary terminals operating in the same spectrum 
band. For the IT-SS model of CRs, the exact operation model 
of the primary network is not important to our study, provided 
that all the secondary terminals satisfy the prescribed IT 
constraints to protect the primary terminals. Without loss of 
generality, we assume that there are K secondary links and J 
primary terminals in each type of the CR networks. 

Consider first the infrastructure-based secondary/CR net- 
work with the S-BS coordinating all the CR transmissions, 
which usually corresponds to one particular cell in a CR 
cellular network. The uplink transmissions from the SUs to 
the S-BS are usually modeled by a multiple-access channel 
(MAC), while the downlink transmissions from the S-BS 
to different SUs are modeled by a broadcast channel (BC). 
For the MAC, the equivalent baseband transmission can be 
represented as 

K 

y = Y, HfeXfc + z (1) 

where y G ([^Mxi denotes the received signal at the S-BS, with 
M denoting the number of antennas at S-BS; Hfe G C^^^^'= 
denotes the channel from the kth SU to S-BS, k ~ 1, ■ ■ ■ , K, 
with Nk denoting the number of antennas at the fcth SU; £ 
([^Nk X 1 (jejjotes the transmitted signal of the fcth SU; and z G 
^A/xi (jgjjotes the noise received at S-BS. We assume that 
Xfe's are independent over fc. 

Similarly, the BC can be represented as 

yfe = Hfx + Zfc, k = l,---,K (2) 

where e c^fcxi denotes the received signal at the fcth 
SU; for convenience, we have used the Hermitian transposed 
uplink channel matrix for the corresponding downlink channel 
matrix, i.e., denotes the channel from the S-BS to the fcth 
SU; X e c^/xi (jgjjotes the transmitted signal from S-BS; 
and Zfc G c^fcXi (jgjjotes the receiver noise of the fcth SU. 
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In the case that x carries information common to all SUs, the 
associated downlink transmission is usually called multicast, 
while if X carries independent information for different SUs, 
it is called unicast. 

Next, consider the ad hoc secondary/CR network, which 
is usually modeled as an interference channel (IC). For 
convenience, we assume that for the fcth secondary link, 
k = I,-- - ,A', the transmitter is denoted as SU-TX^ and 
the receiver is denoted as SU-RX^, although in general a 
secondary terminal can be both a transmitter and a receiver 
The baseband transmission of the IC can be represented as 

K 

Yk = HfefeXfc + ^ HifeXi + Zfc, fc = 1, • • • , A' (3) 

z— 1,27^ A: 

where yk £ 0^*=^^ denotes the received signal at SU-RXfe, 
with Bk denoting the number of receiving antennas; Xfc G 
C^fcxi denotes the transmitted signal of SU-TXfe, with Ak 
denoting the number of transmitting antennas; H^fc G c^fc 
denotes the direct-link channel from SU-TXfe to SU-RXfc, 
while Hife G C^*" ^'^^ denotes the cross-link channel from SU- 
TXi to SU-RXfc, i ^ k; and Zfc G C'^'"^^ denotes the noise at 
SU-RXfc. It is assumed that Xfc's are independent over k. 

Furthermore, we assume that the jth PU, j = 1, • • ■ , J, in 
each type of the CR networks is equipped with Dj antennas, 
Dj > 1. We then use G^j G C^^^^*- to denote the channel 
from the fcth SU to the jth PU in the CR MAC, G C^^^*^ 
to denote the channel from S-BS to the jih PU in the CR BC, 
and Efcj G C^^^'^'' to denote the channel from SU-TXfc to 
the jth PU in the CR IC. Moreover, the receiving terminals 
in the secondary networks may experience interference from 
active primary transmitters. For simplicity, we assume that 
such interference is treated as additional noise at the secondary 
receivers, and the total noise at each secondary receiving 
terminal is distributed as a CSCG random vector with zero 
mean and the identity covariance matrix. 

Note that the (spatial) channels defined in the above CR 
network models are assumed constant for a fixed transmit 
dimension such as one time-block in a time-division-multiple- 
access (TDMA) system or one frequency-bin in an orthogonal- 
frequency-division-multiplexing (OFDM) system. In a wire- 
less environment, these channels usually change over time 
and/or frequency dimensions as governed by an underlying 
joint stochastic process. As such, DRA becomes relevant to 
schedule SUs into different transmit dimensions based on 
their CSI. In general, the secondary transmitting terminals 
need to satisfy two types of power constraints for DRA: 
One is due to their own transmit power budgets; the other 
is to limit their resulting interference level at each PU to be 
below a prescribed threshold. These constraints can be applied 
over each fixed dimension as peak power constraints, or over 
multiple dimensions as average power constraints. Without 
loss of generality, we consider DRA for the secondary network 
over L transmit dimensions with different channel realizations, 
with L > 1. In total, four different types of power constraints 
can be defined for the secondary network. By taking the CR 
MAC as an example (similarly as for the CR BC/IC), we have 



. Peak transmit power constraint (PTPC): 

Tr(Sfe[/]) < Pfc (4) 

where Sfc [I] denotes the transmit covariance matrix for the 
Zth transmit dimension of the fcth SV, I G {1, ■ ■ ■ , L},k € 
{1, • • • ,K}; Pfe denotes the fcth SU's peak power con- 
straint that applies to each of the L transmit dimensions. 
. Average transmit power constraint (ATPC): 

1 ^ 

-^Tr(SfcH)<Pfc (5) 

1=1 

where Pk denotes the fcth SU's average transmit power 
constraint over the L transmit dimensions. 

• Peak interference power constraint (PIPC): 

K 

5^Tr(Gfc,[Z]Sfc[;]Gg.W) <r, (6) 

k=l 

where Gkj [I] denotes the realization of channel Gkj for 
a given I; and Tj denotes the peak interference power 
constraint for protecting the jth PU, j G {I,-- - , J}, 
which limits the total interference power caused by all 
the K SUs across all the receiving antennas of the jth 
PU, for each of the L transmit dimensions. 

• Average interference power constraint (AIPC): 

L K 

L^T.^' {GkAl]Skmg[l]) < f, (7) 

1=1 k=l 

where Tj denotes the average interference power con- 
straint for the jth PU to limit the total interference power 
from the K SUs, which is averaged over the L transmit 
dimensions. 

Note that DRA for traditional (non-CR) wireless networks 
under PTPC and/or ATPC has been thoroughly studied in the 
literature (TTl . while the study of DRA subject to PIPC and/or 
AIPC as well as their various combinations with PTPC/ ATPC 
is unique to CR networks and is relatively new. In order to gain 
more insights into the optimal DRA designs for CR networks, 
we will first study the case of a single transmit dimension 
(L = 1) with PTPCs and PIPCs by focusing on the spatial- 
domain transmit optimization for multi-antenna CRs in Section 
Hin and then study the general case of L > 1 for joint space- 
time-frequency DRA in CR networks under ATPCs and AIPCs 
in Section |IV] 

III. Cognitive Beamforming Optimization 

In this section, we consider the case of L = 1, where the 
DRA for CR networks reduces to the spatial-domain transmit 
optimization under PTPCs and PIPCs to maximize the CR 
network throughput. We term this practice as cognitive beam- 
forming. In order to investigate the fundamental performance 
limits of cognitive beamforming, we study the optimal designs 
with the availability of perfect knowledge on all the channels 
in the SU networks, and those from all the secondary transmit 
terminals to PUs. For convenience, we drop the dimension 
index I for the rest of this section given L = 1. 
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First, it is worth noting that the PIPC given in (|6]) can be 
unified with the PTPC given in ^ into a form of generalized 
linear transmit covariance constraint (GLTCC): 

K 

5^Tr(W,S0<u. (8) 

i=l 

where W^'s and if are constants. For example, with each PIPC 
given in (|6]), = Gf^Gy , Vi, and w ~ Tj, while for each 
PTPC given in (HJi, Wj = 1 \f i ~ k and otherwise, with 
w ^ Pk- Previous studies on transmit optimization for multi- 
antenna or multiple-input multiple-output (MIMO) systems 
have mostly adopted some special forms of GLTCC such as the 
user individual power constraints and sum-power constraint. 
However, it remains unclear whether such existing solutions 
are applicable to the general form of GLTCC, which is crucial 
to the problem of CR MIMO transmit optimization with the 
newly added PIPCs. In the following, we provide an overview 
of the state-of-art solutions for this problem under different 
CR network models, while the developed solutions also apply 
to the case with the general form of GLTCCs as in ([8]). 
From a convex optimization perspective, we next divide our 
discussions into two parts, which deal with the cases of convex 
and non-convex problem formulations, respectively. 

A. Convex Problem Formulation 

First, consider the case where the associated optimization 
problem in a traditional MIMO system without PIPC is convex. 
In such cases, since the extra PIPCs are linear over the SU 
transmit covariance matrices, the resulting transmit covariance 
optimization problem for CR systems remains convex; and 
thus, it can be efficiently solved by standard convex optimiza- 
tion techniques. 

CR Point-to-Point MIMO Channel: We elaborate this 
case by first considering the CR point-to-point MIMO channel, 
which can be treated as the special case with only one active 
SU link in the MAC, BC, or IC based CR network. Without 
loss of generality, we will use the notations developed for 
the CR MAC with if = 1 in the following discussions. 
Specifically, the optimal transmit covariance to achieve the CR 
point-to-point MIMO channel capacity under both the PTPC 
and PIPCs can be obtained from the following problem |fT3]| : 

Max. log|l + HSH^| (PI) 

s. t. Tr(S) < P 

Tr(G,SGf) <r„ i = l,-.. ,J 
S h 

where for conciseness we have removed the SU index k in the 
symbol notations since K = \, while S ^ means that S is 
a positive semi-definite matrix lfT4ll . 

We see that (PI) is a convex optimization problem since 
its objective function is concave over S and its constraints 
define a convex set over S. Thus, (PI) can be efficiently 
solved by, e.g., the interior point method lfT4l . In the special 
case of CR multiple-input single-output (MISO) channel, i.e., 
H degrades to a row-vector denoted by h G C^^^, it 
can be shown by exploiting the Karush-Kuhn-Tucker (KKT) 



optimality conditions of (PI) that transmit beamforming is 
capacity optimal, i.e., Rank(S) = 1 |fT3l|. Thus, without loss 
of generality we could write S = vv^, where v G C^^^ 
denotes the precoding vector. Accordingly, (PI) for the special 
case of MISO CR channel is simpHfied as (Pl-S) HI: 

Max. ||hv|| 

V 

s. t. ||vf <P 

l|G,v||2<rj, - ,j, 

which is non-convex due to the non-concavity of its objective 
function. However, by observing the fact that if v is the 
solution of (Pl-S), so is e-'^v for any arbitrary 0, we thus 
assume without loss of generality that hv is a real number and 
modify (Pl-S) by rewriting its objective function as Rc(hv) 
and adding an additional linear constraint Im(hv) = 0. 
Thereby, (Pl-S) can be converted into a second-order cone 
programming (SOCP) 1141 problem, which is convex and thus 
can be efficiently solved by available convex optimization 
softwares ifTsll . Alternatively, (Pl-S) can be shown equiva- 
lent to its Lagrange dual problem 1131 . which is a convex 
semi-definite programming (SDP) lfT4l problem and is thus 
efficiently solvable 031 . For (Pl-S) in the case of one single- 
antenna PU, a closed-form solution for the optimal precoding 
vector V was derived in ifTsl via a geometric approach. 

In order to reveal the structure of the optimal S for (PI), 
we consider its Lagrange dual problem defined as (Pl-D): 

Mill, diri) 

where 77 = [?7o7 • • • , "hA denotes a vector of dual variables 
for (PI) with 770 associated with the PTPC, and rjj associated 
with the J th PIPC, j = 1 , • ■ • , J, while we have the dual 
function defined as 

d{r]) ^ max log |l + HSH^| - 7?o(Tr(S) - P) 
J 

-5]ry,(Tr(G,SGf)-r,). (9) 

Since (PI) is convex with Slater's condition satisfied lfT4l . 
the duality gap between the optimal values of (PI) and (PI- 
D) is zero, i.e., (PI) can be solved equivalently as (Pl-D). 
Accordingly, an iterative algorithm can be developed to solve 
(Pl-D) by alternating between solving d{ri) for a given r] and 
updating r] to minimize d{ri). At each iteration, r] can be 
updated by a subgradient-based method such as the ellipsoid 
method [16|, according to the subgradients of d{ri), which can 
be shown equal to P-Tr(S*) and Tj-Tr {GjS*Gf) for tjq 
and rjj ,j^0, respectively, where S* denotes the optimal S to 
obtain d{ri) for a given rj. From (|9]), it follows that S* is the 
solution of the following equivalent problem (by discarding 
irrelevant constant terms): 

max logll + HSH^I -Tr(TS) (10) 
s^o ' I \ / 

where T = 770! + J2j=i Vji^j^^j) is ^ constant matrix for 
a given rj. In order to solve Problem ( fTOl i. we introduce an 
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auxiHary variable: S = tI/^STI/^. Problem ^ is then re- 
expressed in terms of S as 



max log 
s^o 



I + HT"1/2st-i/2h^ 



Tr(S). 



(11) 



The above problem can be shown equivalent to the standard 
point-to-point MIMO channel capacity optimization problem 
subject to a single sum-power constraint flTll . and its solution 
can be expressed as S* = VSV^, where V is obtained from 
the singular-value decomposition (SVD) given as follows: 
jjT-i/2 = U0V^, with = Diag( [6*1,..., 6It]) and 
T = min{M,N), while S = Diag([cri, . . . , ctt]) follows 
the standard water- filling solution ifTTl : ai = (l/ln2 — 
1/6*2)+, i = 1,...,T, with (■)+ = max(0, •)■ Thus, the 
solution of Problem ( fTOl i for a given t] can be expressed as 

Next, we present a heuristic method for solving (PI), which 
leads to a suboptimal solution in general and could serve 
as a benchmark to evaluate the effectiveness of the above 
two approaches based on convex optimization. To gain some 
intuitions for this method, we first take a look at two special 
cases of (PI). For the first case, supposing that all the PlPCs 
are inactive (e.g., by setting Tj = oo,Vj) and thus can be 
removed, (PI) reduces to the standard MIMO channel capacity 
optimization problem under the PTPC only, for which the 
optimal solution of S is known to be derivable from the SVD 
of H |17l|. For the second case, assuming that Tj = 0,Vj', 
the solution for (PI) is then obtained by the "zero-forcing 
(ZF)" algorithm ITSl . which first projects H into the space 
orthogonal to all Gj's, and then designs the optimal S based 
on the SVD of the projected channel. Note that the (non- 
trivial) ZF-based solution exists only when N > X]/=i ^j- 
From the above two special cases, we observe that as Fj's 
decrease, the optimal S should evolve along with a sequence 
of subspaces of H with decreasing dimensions as a result of 
keeping certain orthogonality to Gj 's, which motivates a new 
design method for cognitive beamforming, named as partial 
channel projection ifTsl . Specifically, let Gj = Gj/Tj,yj. 
Then, define G = [Gf , • • • ,G^f. Denote the SVD of G 
as G = UgAgVq. Without loss of generality, assume that 
the singular values in Aq arranged in a decreasing order. 
Then, we propose a generalized channel projection operation: 



H 



HIV 



V 



(b) 



H 



(12) 



where Vq-* consists of the first b columns of Yq corre- 
sponding to the b largest singular values in Ag, 1 < 6 < 
min(A^ — 1, J2 j=i ^j)- Note that b could also take a zero value 
for which Yq = 0. Now, we are ready to present the transmit 
covariance matrix for the partial projection method in the form 
of its eigenvalue decomposition (EVD) as S = Vj^Sj^V^, 
where Vj^ is obtained from the SVD of the projected channel 
Hj^, i.e., = U_lAj^V^. By substituting this new form 
of S into (PI), it can be shown that the problem reduces to 
maximizing the sum-rate of a set of parallel channels (with 
channel gains given by A^) over their power allocation 'E± 
subject to (J + 1) linear power constraints, for which the 
optimal power allocation can be obtained by a generalized 



"water-filling" algorithm lfT3l . Note that the partial channel 
projection works for any values of N and Dj's. 

In Fig. m we plot the achievable rate of a CR MIMO channel 
under the PTPC and PlPCs with the optimal transmit covari- 
ance solution for (PI) via the convex optimization approach, 
against those with suboptimal covariance solutions via the 
partial channel projection method with different values of b. 
The system parameters are given as follows: M = N = 4, 
J 2, Di = 152 1, and Fi = Fa = 0.1. The 
SU achievable rate is plotted vs. the SU PTPC, P. It is 
observed that the optimal covariance solution obtained via 
the convex optimization approach yields notable rate gains 
over suboptimal solutions via the heuristic method, for which 
the optimal value of b (the number of SU-to-PU channel 
dimensions to be nulled) to maximize the SU achievable rate 
increases with the SU PTPC. 

CR MIMO-MAC: The solutions proposed for the CR 
point-to-point MIMO channel shed insights to transmit opti- 
mization for the CR MIMO-MAC defined in O with K > 1. 
Assume that in the CR MIMO-MAC, the optimal multiuser 
detection is deployed at the S-BS to successively decode 
different SU messages from the received sum-signal. We then 
consider the problem for jointly optimizing SU transmit co- 
variance matrices to maximize their weighted sum-rate subject 
to individual PTPCs and joint PlPCs. This problem is referred 
to as weighted sum-rate maximization (WSRMax). Without 
loss of generality, we assume that the given user rate weights 
satisfy that /ii > /i2 > • ■ • > l^K > 0; thus, the optimal 
decoding order of users at the S-BS to maximize the weighted 
sum-rate is in accordance with the reverse user index lfT9l . 
Accordingly, the WSRMax for the CR MIMO-MAC can be 
expressed as 

K I + ^,tiH.S,Hf 
Max. V ^lh log (P2) 

s. t. Tr(Sfe) <Pfc, k^l,...,K 

K 

^Tr(Gfc,SfcGg-) <F„ j^l,---,J 

k=l 

Sfc ^ 0, k^l,--- ,K. 
Reordering terms in the objective function of (P2) yields 



K-l 



k=l 



I + ^ H,S,Hf 



IJ-K log 



K 



I + ^H,S,Hf 



(13) 



From the above new form of the objective function, it can 
be verified that (P2) is a convex optimization problem over 
Sfc's. Thus, similarly as for (PI), (P2) can be solved by an 
interior-point-method based algorithm or an iterative algorithm 
via solving the equivalent Lagrange dual problem, for which 
the details are omitted here for brevity. 

It is noted that (P2) is for the case with the optimal non- 
hnear multiuser decoder at the S-BS, while in practice the 
low-complexity linear decoder is usually more preferable. The 
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- Optimal Solution 

- Suboptimal Solution {b=0} 

- Suboptimal Solution {b=1 

- Suboptimal Solution {b=2} 
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Fig. 2. Comparison of the achievable rates for the CR MliVIO channel under 
the PTPC and PIPCs with the optimal transmit covariance solution for (PI) 
via the the convex optimization approach vs. suboptimal covariance solutions 
via the partial channel projection method with different values of b. 

use of linear instead of non-linear decoder at the receiver 
will change the user achievable rates for the CR MIMO- 
MAC, thus resulting in new problem formulations for transmit 
optimization. For example, in 1201 . the authors have considered 
the CR SIMO-MAC (single-antenna for each SU transmitter) 
with a Unear decoder at the receiver, where the power alloca- 
tion across the SUs is optimized to maximize their signal-to- 
interference-plus-noise ratios (SINRs) at the receiver subject 
to both transmit and interference power constraints. 



B. Non-Convex Problem Formulation 

Next, we consider the case where the optimization problems 
in the associated traditional (non-CR) MIMO systems are non- 
convex. It thus becomes more challenging whether these non- 
convex problems with the addition of convex PIPCs in the 
corresponding CR MIMO systems can be efficiently solvable. 
In the following, we present some promising approaches to 
solve these problems for the CR MIMO-BC and MIMO-IC. 

CR MIMO-BC: First, consider the CR MIMO-BC defined 
in ^ under both the PTPC at the S-BS and J PIPCs each 
for one of the J PUs, which can be similarly defined as for 
the MAC case in (HJi and (|6]l, respectively. We focus on the 
unicast downlink transmission for the CR BC, while for the 
case of multicast, the interested readers may refer to 121)1 . For 
the purpose of exposition, we consider two commonly adopted 
design criteria for the traditional multi-antenna Gaussian BC 
in the literature: One is for the MIMO-BC deploying the 
non-linear "dirty -paper-coding (DPC)" at the transmitter l22ll . 
which maximizes the weighted sum-rate of all the users 
(i.e., the WSRMax problem); the other is for the MISO-BC 
(single-antenna for each SU receiver) deploying only linear 
encoding at the transmitter, which maximizes the minimum 
SINR among all the users, referred to as "SINR balancing". 

Specifically, the WSRMax problem for the CR MIMO-BC 
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can be formulated as: 
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Max ^MUog 
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s. t. ^Tr(Qfe)<F 

fc=i 

Tr(^F, (j2Clk^ Ff^ <r„ j = l,... ,J 

Qfc^O, k = l,--- ,K 

where e (^mxm (jgnotes the covariance matrix for 
the transmitted signal of S-BS intended for the kth SU, 
k = 1, ■ • • ,K; Hk's. are the given user rate weights; and P 
denotes the transmit power constraint for the S-BS. Without 
loss of generality, we assume that /ii > /i2 > • • ■ > I^k > 0; 
thus, in (P3) the optimal encoding order of users for DPC 
to maximize the weighted sum-rate is in accordance with 
the user index l22l . Note that (P3) is non-convex with or 
without the PIPCs due to the fact that the objective function 
is non-concave over Q^'s for iiT > 2. As a result, unlike 
(PI) for the point-to-point CR channel, the standard Lagrange 
duality method cannot be applied for this problem. For (P3) 
in the case without the PIPCs, a so-called "BC-MAC duality" 
relationship was proposed in f23l to transform the non-convex 
MIMO-BC problem into an equivalent convex MIMO-MAC 
problem, which is solvable by efficient convex optimization 
techniques such as the interior point method. In l24l . another 
form of BC-MAC duality, the so-called "mini-max duality" 
was explored to solve the MIMO-BC problem under a special 
case of GLTCC: the per-antenna transmit power constraint. 
However, these existing forms of BC-MAC duality are yet 
unable to handle the case with arbitrary numbers of GUTCCs, 
which is the case for (P3) with both the PTPC and PIPCs. 

In l25]| . a general method was proposed to solve various 
MIMO-BC optimization problems under multiple GLTCCs, 
thus including the CR MIMO-BC WSRMax problem given in 
(P3). For this method, the first step is to combine all (J + 1) 
power constraints in (P3) into a single GLTCC as shown in 
the following optimization problem: 



K 

Max. > iih loe 



fe=i 



I + Hf (EL+iQOH 



s. t. 



Qfc ^ 0, fc = I,-- - ,A' 



(14) 



where A = AoI+E/=i AjFf Fj, and Q = AoP+E/=i >^3^o 
with Ao, Ai, • • • , Aj being non-negative constants. For a given 
set of Ai's, i = 0, ■ • • , J, let the optimal value of the 
above problem be denoted by F(Ao,Ai,-- - ,Aj). Clearly, 
i^(Ao,Ai,-- - ,Aj) is an upper bound on the optimal value 
of (P3) since any feasible solutions for (P3) must satisfy 
the constraints of Problem ( fT4b for a given set of A/s. 
Interestingly, it can be shown that the optimal value of (P3) 
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Fig. 3. Generalized MIMO MAC-BC Duality: (a) Primal MIMO-BC channel 
with downlink channels and receiver noise vectors zj. ~ CA^(0, 1), 

k = 1,...,K, and a GLTCC: Tr{AY,k=i Qk) < Q\ (b) Dual MIMO- 
MAC with uplink channels Hfc,A: = 1,...,K and receiver noise vector 
z ~ CA^(0,A), and a sum-power constraint: I]^iTr(Sfe) < Q. The 
MIMO-BC and dual MIMO-MAC have the same achievable rate region l25l . 

is equal to the minimum value of function -F(Ao, Ai, • ■ • , Aj) 
over all non-negative A^'s ||25l . Therefore, (P3) can be re- 
solved by iteratively solving Problem (fl4] i for a given set 
of Ai's and updating A/s towards their optimal values to 
minimize function F(Ao, Ai, ■ • ■ , Aj). Specifically, A/s can be 
updated via the ellipsoid method according to the subgradients 
of F(Ao,Ai,-- - ,Aj), which can be shown 1251 equal to 
P - Eti Tr(Q^) and - Tr(F,(E[ii Q^)Ff ) for A„ 
and Xj (j 0), respectively, where Q^'s are the solution of 
Problem (fT4] i for the given Afe's. 

Furthermore, Problem (fl4] i with a given set of A^'s can be 
solved by applying the generalized BC-MAC duality proposed 
in ||251 . which extends the existing forms of BC-MAC duality 
||23l . Il24l to transform the MIMO-BC problem subject to 
a single GLTCC as in Problem (fl4] i to an auxiliary (dual) 
MIMO-MAC problem subject to a corresponding sum-power 
constraint. Specifically, it is shown in [25 1 that the MIMO- 
BC as in Problem (O and the dual MIMO-MAC, as depicted 
in Fig. |3] have the same achievable rate region. Accordingly, 
the optimal objective value (weighted sum-rate) of Problem 
(0 for the primal MIMO-BC can be obtained as that of the 
following equivalent problem for the dual MIMO-MAC: 



Max. 
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fc=i 
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s. t. ^Tr(Sfc)<g 

k=l 

Sfe ^ 0, fc - 1, ■ • 



(15) 



Similar to (P2), the above problem is a WSRMax problem for 
the MIMO-MAC subject to a single sum-power constraint, 
which is convex and thus can be efficiently solvable by, e.g., 
the interior point method. After solving Problem ( fTsT i, the 
optimal user transmit covariance solutions for the MIMO- 
MAC, S^'s, can be transformed to the corresponding ones 
for the original MIMO-BC, Q^'s, via a MAC-BC covariance 
transformation algorithm given in 1251 . Furthermore, it is 



worth noting that with K = 1, the above method can be shown 
equivalent to that developed for (PI) in the CR point-to-point 
MIMO channel case based on the Lagrange duality. 

Consider next the SINR balancing problem for the CR 
MISO-BC, which can be expressed as: 



Max. a 

a,vi,--- .V K 
S. t. 



|hfv,f 



(P4) 

> a, k = 1, - ■ • ,K 



k=l 
K 



fc=i 

where a denotes an achievable SINR for all the SUs; G 
(j-<Mxi (jgjjotes the precoding vector for the transmitted signal 
of S-BS intended for the fcth SU; and represents for 
the MISO-BC case. Similarly as for (Pl-S), by treating Vfc 
on the left-hand side (LHS) of each SINR constraint in (P4) 
as a positive real number ll26l . it can be shown that (P4) for a 
given a is equivalent to a SOCP feasibility problem and thus 
efficiently solvable ifTSl . For a given a, if the associated SOCP 
problem is feasible, we know that the optimal solution of (P4) 
for a, denoted by a*, must satisfy a* > a; otherwise, a* < a. 
Based on this fact, a* can be found by a simple bisection 
search lT4l : with a*, the corresponding optimal solution for 
Vfc's in (P4) can also be obtained. The above technique has 
also been applied in ll27l for (P4) without the PIPCs. 

The SINR balancing problem for the conventional MISO- 
BC without the PIPCs has also been studied in [281, where an 
algorithm was proposed using the virtual uplink formulation 
and a fixed-point iteration. However, this algorithm cannot be 
extended directly to deal with multiple PIPCs for the case of 
CR MISO-BC. Similarly as for the previous discussions on the 
WSRMax problem for the CR MIMO-BC where a generalized 
MIMO MAC-BC duality holds, a counterpart beamforming 
duahty also holds for the MISO-BC and SIMO-MAC l25l . 
With this duality result, the SINR balancing problem (P4) 
for the CR MISO-BC can be converted into an equivalent 
problem for the dual SIMO-MAC, where the efficient iterative 
algorithm in ll28l can be directly applied. The interested 
readers may refer to ll25l for the details of this method. 

CR MIMO-IC: Second, consider the CR MIMO-IC given 
in (O), subject to both the PTPCs for the K SU-TXs and the 
PIPCs for the J PUs, which can be similarly defined as for the 
MAC case in (|4]i and (|6]l, respectively. From an information- 
theoretic perspective, the capacity region for the Gaussian IC 
under PTPCs, which consists of all the simultaneously achiev- 
able rates of all the users, still remains unknown in general 
even for the case of K = 2 and Ai. ^ Bi, = I. k = 1,2 l29l. A 
pragmatic approach that leads to suboptimal achievable rates 
in the Gaussian IC is to restrict the system to operate in a 
decentralized manner, i.e., allowing only single-user encoding 
and decoding by treating the co-channel interferences from the 
other users as additional Gaussian noises. For this approach, 
transmit optimization for the CR MIMO-IC reduces to finding 
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a set of optimal transmit covariance matrices for the K SU 
links, denoted by 6 C-^^ x — i . . . , /v , to maximize 
the secondary network throughput under both the PTPCs and 
PIPCs. More specifically, the WSRMax problem for the CR 
MIMO-IC can be expressed as: 



K 



k=l 



Max ^ Mfc log I + (I + E H.fcR.H^)-iHfc,.RfeHf, 

(P5) 



s. t. Tr(Rfc) <Pfc, = - ,if 
if 

ETr(Efc,RfcEg) <r,, j = l,-- - ,J 

k=l 

Rfe ^0, fc = l,-- - ,if 

where /i^'s are the given non-negative user rate weights. We 
see that (P5) is non-convex with or without the PIPCs due to 
the fact that the objective function is non-concave over R^'s 
for K > 1. As a result, there are no efficient algorithms yet 
to obtain the globally optimal solution for this problem. For 
the same problem setup, there have been recent progresses on 
characterizing the maximum achievable "degrees of freedom 
(DoF)" for the user sum-rate (i.e., /Zfc = l,V/c) ll30l . 

Next, we discuss some feasible solutions for (P5). First, it 
is worth noting that for (P5) in the case without the PIPCs, a 
commonly adopted suboptimal approach is to iteratively opti- 
mize each user's transmit covariance subject to its individual 
PTPC with the transmit covariances of all the other users 
fixed. This decentralized approach has been first proposed 
in II3TI . 1321 to obtain some local optimal points for (P5) 
with the PTPCs only, where they differ in that the one in 
II3TI maximizes the user individual rate at each iteration, 
while the one in |32l maximizes the user weighted sum-rate. 
It is also noted that a parallel line of works with similar 
iterative user optimizations has been pursued in the single- 
antenna but multi-carrier based interference channels such as 
the wired discrete-multi-tone (DMT) based digital subscriber 
line (DSL) network lf33l , and the wireless OFDM based ad 
hoc network 1341 . One important question to answer for such 
iterative algorithms is under what conditions the algorithm will 
guarantee to converge to a local optimal point. This problem 
has been addressed in the contexts of both multi-carrier and 
multi-antenna based interference channels in, e.g., ||35l , l36l , 
via game-theoretic approaches. 

However, the above iterative approach cannot be applied 
directly to solve (P5) with both the PIPCs and PTPCs, since 
each PIPC involves all the user transmit covariances and is 
thus not separable over the SUs. Thus, a feasible approach 
for (P5) is to decompose each of the J PIPCs into a set of 
interference-power constraints over the K SU-TXs, i.e., for 
the jth PIPC, j e {I,-- - , J}, 



Tr fEfc,Rj.EE) < r 



kj) 



fc = 1, 



(16) 



where F^*^'' is a constant, and all F^'^^'s, k ~ 1 

satisfy J^k ^j''^ — ^^^^ that the jth PIPC is guaranteed. 
Then, the iterative algorithm works here, where each SU link 
independently optimizes R^ to maximize its achievable rate 



under its PTPC and J interference-power constraints given by 
(fTSI l. with all other R^'s, i ^ k, fixed. It is observed that 
the resulting problem is in the same form of our previously 
studied (PI) for the CR point-point MIMO channel; thus, 
similar techniques developed for (PI) can be applied. Note 
that a suboptimal method for this problem in the same spirit 
of the partial channel projection method to reduce the design 
complexity for each SU transmit covariance matrix has also 
been proposed in ll37l . Moreover, it is noted that F^'^^'s, 
j ~ 1, • ■ • , J, k = 1, ■ • • ,K, can be searched over the SUs 
to further improve their weighted sum-rate. 

Alternatively, assuming that a centralized optimization is 
feasible with the global knowledge of all the channels in the 
SU network, as well as those from different SU-TXs to all 
PUs, another heuristic algorithm for (P5) was proposed in l38l . 
By rewriting the SU transmit covariance matrices into their 
equivalent precoding vectors and power allocation vectors, 
this algorithm iteratively updates the SU transmit precoding 
vectors (based on the "network duality" l39l ) or the power 
allocation vectors (by solving geometric programming (GP) 
problems ||40J ), with the others being fixed. 

It is worth pointing out that there are other problem formu- 
lations different from (P5) to address the transmit optimization 
for the CR MIMO-IC. In BTI . a new criterion was proposed 
to design the SU link transmission in a CR MISO-IC via an 
alternative decentralized approach, where each SU-TX inde- 
pendently designs its transmit precoding vector to maximize 
the ratio between the received signal power at the desired SU- 
RX and the resulted total interference power at all the PUs, 
in order to regulate the interference powers at PUs. Moreover, 
the above discussions are all based on the assumption that 
each SU-RX treats the interferences from all the other SU 
links as additional noises, which is of practical interest since 
it simplifies the receiver design for each SU link. However, 
due to independent cross-link channels between SU terminals, 
it may be possible that a SU-RX could occasionally observe 
"strong" interference signals from some co-existing SU-TXs 
and thus be able to decode their messages via multiuser detec- 
tion techniques and then cancel the associated interferences. 
With such "opportunistic" multiuser detection at each SU-RX, 
the achievable rate of each SU link becomes a function of 
not only its own transmit covariance, but also those of the 
other SUs as well as their instantaneous transmit rates. Thus, 
the corresponding transmit optimization for the CR MIMO-IC 
leads to new and more challenging problem formulations than 
(P5); the interested readers may refer to 1421 . l43l . 

IV. Joint Space-Time-Frequency DRA Optimization 

In the previous section, we have studied DRA for different 
CR networks at a single transmit dimension in time/frequency, 
by focusing on spatial-domain transmit optimization under the 
peak transmit and interference power constraints (PTPC and 
PIPC). In this section, we bring the additional time and/or 
frequency dimensions into the DRA problem formulations, 
by applying the average transmit and interference power con- 
straints (ATPC and AIPC) in CR networks. Consider the DRA 
over L time/freqeuncy dimensions, for which all the required 
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channel knowledge is assumed to be known. Taking the CR 
MAC as an example (similar arguments can be developed for 
the CR BC/IC), under both the ATPCs and AIPCs given in ^ 
and (|7]i, respectively, a generic problem formulation for DRA 
optimization can be formulated as: 



Max. C({SfcW}) 

Sfc[i]^0,Vfc,/ 



(P6) 



s. t. (El,© 



where {Sfe[Z]} denotes the set of Sfc[/]'s, k = I,-- - ,K, 
and I = I,-- - ,L, while C(-) is an arbitrary utility func- 
tion to measure the CR network performance. We as- 
sume that C(-) is separable over Z's, i.e., C({Sfc[^]}) = 
T Ef=i Ui{Si[l], • • • , Sk[1]) with Ui{-ys denoting individual 
utility functions. Since both the ATPC and AIPC involve L 
transmit covariance matrices, the Lagrange dual decomposi- 
tion (see, e.g., a tutorial paper ll44ll ) is a general method to deal 
with this type of average constraints for optimization over a 
number of parallel dimensions, which is explained as follows. 
By introducing a set of dual variables, z/fc's, each for one of 
the K ATPCs, and Sk's, each for one of the J AIPCs, the 
Lagrange dual problem of (P6) can be written as (P6-D): 



Min. 



with V = [i^i, • ■ • , i^k], 6 = [6i, ■ ■ ■ ,6j], and the dual 
function 



d{u,S) ^ max C({SfeH}) - ^ ^.,(1 ^ Tr(Sfc[; 

k=l 1=1 

.7 L K 



j = l 1=1 k=l 



(17) 



Since the dual problem (P6-D) is convex regardless of the 
convexity of the primal problem (P6) il4] . (P6-D) can be 
efficiently solved by the ellipsoid method according to the 
subgradients of the dual function d{v, S), similarly as in our 
previous discussions, provided that the maximization problem 
in ([TtI i is solvable for any given set of 1/ and 6. It is interesting 
to observe that this maximization problem can be decomposed 
into L parallel subproblems each for one of the L dimensions, 
and all of these subproblems have the same structure and 
are thus solvable by the same algorithm, a practice known 
as "dual decomposition". Without loss of generality, we drop 
the dimension index I and express each subproblem as 



K 



max U(Si, 
Skhoyk 



,SK)-ETr(BfcK,^)Sfc) 



(18) 



fc=i 



where Bk{i^k,S) = Vk^ + SjLil'^jG^Gfcj) is a constant 
matrix for the given Vk and S, k = 1 ■ ■ ■ ,K. 

We then discuss the following two cases. For the first case, 
consider that [/;(•) is a concave function over Sfc[/]'s, 
(e.g., the point-to-point CR channel capacity in (PI), or the 
weighted sum-rate for the CR MIMO-MAC in (P2)). Then, 
(P6) is convex and thus the duality gap between the optimal 
values of (P6) and (P6-D) is zero, i.e., (P6) and (P6-D) are 



equivalent problems. Furthermore, each subproblem in (fTsT l 
is also convex. Thus, the dual decomposition method solves 
(P6) via its dual problem (P6-D), which is decomposable into 
L convex subproblems. For the second case, as a counter- 
part, consider that [/;(•) is non-concave over Sfe[Z]'s (e.g., 
the weighted sum-rate for the CR MIMO-BC/MIMO-IC in 
(P3)/(P5)). As a result, (P6) is non-convex and the duality gap 
between (P6) and (P6-D) may not be zero. Furthermore, the 
subproblem (fTSl l is also non-convex. For this case, even when 
the optimal solutions of the L subproblems are obtainable, the 
optimal value of (P6-D) in general only serves as an upper 
bound on that of (P6). However, in [45] it is pointed out that 
if a set of so-called "time-sharing" conditions are satisfied by 
a non-convex optimization problem, the duality gap for this 
problem and its dual problem is indeed zero. Furthermore, for 
the class of DRA problems in the form of (P6), the associated 
time-sharing conditions are usually satisfied asymptotically 
as L — > 00 under some cautious considerations on the 
continuity of channel distributions P6| . Therefore, the dual 
decomposition method could still be applied to solve (P6) in 
the non-convex case for sufficiently large values of L, provided 
that the optimal solutions for the subproblems in ( fTsT i are 
obtainable (e.g., a variation of (P3) for the CR MIMO-BC). 
However, with finite values of L, how to efficiently solve (P6) 
in the case of non-concave objective functions is still open. 

With the above discussions on the general approaches to 
design joint space-time-frequency DRA for CR networks, we 
next present some examples of unique interests to CR systems. 



A. TDMA/FDMA Constrained DRA: When Is It Optimal? 

Time-/frequency-division multiple-access (TDMA/FDMA), 
which schedules only one user for transmission at each 
time/frequency dimension, is usually preferable in practice due 
to their implementation ease. For the TDMA/FDMA based 
CR MAC (similar arguments hold for the CR BC/IC), the 
optimal DRA over L transmit dimensions to maximize the 
sum-capacity of the SUs can be formulated as (P6) with 
properly chosen functions for [/;(-)'s, where for any given /, 
[/;(•) is expressed as (l is dropped for conciseness) 



[/(Si,-.. ,Sk) = 



log 




I + HfcSfcH 



S, ^ 0,Vi # k 
otherwise. 

(19) 



Note that U{-) defined above implies the TDMA/FDMA 
constraint, i.e., only scheduling one user for transmission at 
a given dimension with a positive contribution to the sum- 
capacity. However, it can be shown that [/(.) is non-concave 
over Sfe's in this case and as a result, the corresponding (P6) 
is non-convex. Nevertheless, according to our previous discus- 
sions, since the time-sharing conditions hold approximately 
when L — !• 00, the dual decomposition method can be applied 
to solve (P6) for this case with very large values of L, where 
the optimal solution of the associated subproblem at each 
dimension given in (fTsl l can be obtained by finding the SU 
(selected for transmission) with the largest objective value of 
the following problem (which is of the same form as Problem 
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( [Tol l and thus solvable in a similar way): 

max logll + HfcSfeHf I -Tr(BfcK-,<5)Sfc). (20) 

An important question to investigate for TDMA/FDMA 
based DRA is how much the performance is degraded as 
compared with the optimal DRA that allows more than one 
users to transmit at a given dimension. From an information- 
theoretic viewpoint, it is thus pertinent to investigate the 
conditions for the optimality of TDMA/FDMA, i.e., when 
they are optimal to achieve the system sum-capacity. For the 
traditional single-antenna fading MAC under the user ATPCs 
over time, it has been shown in BtI that TDMA is optimal 
for achieving the ergodic/long-term sum-capacity. This result 
has been shown to hold for the fading CR MAC and CR BC 
under both the ATPCs and AIPCs in ||48i|. where by exploiting 
the KKT optimality conditions of the associated optimiza- 
tion problems, the optimality conditions for TDMA in other 
cases of combined peak/average transmit/interference power 
constraints have been characterized. For the traditional single- 
antenna IC with interference treated as noise, the optimality 
of TDMA/FDMA for the sum capacity has been investigated 
under the ATPCs in 1491 . ifSOl . It would be interesting to extend 
these results to the case of CR IC under the additional PIPCs 
and/or AIPCs. 

B. Peak vs. Average Interference Power Constraints: A New 
Interference Diversity 

From a SU's perspective, it is obvious that the ATPC/AIPC 
is more flexible than the PTPC/PIPC for DRA under the same 
power threshold and thus results in a larger SU link capac- 
ity. However, from a PU's perspective, it remains unknown 
whether the AIPC or PIPC causes more PU link performance 
degradation. Intuitively speaking, the PIPC should be more 
favorable than the AIPC since the former limits the interfer- 
ence power at the PU to be below certain threshold at each 
time/frequency dimension, while the latter results in varia- 
tions of interference power levels over different dimensions 
although their average level is kept below the same threshold 
as that for the PIPC. 

Somehow surprisingly, in |5T] it is shown that for the 
single-antenna PU fading channel subject to the interference 
from a SU transmitter, the AIPC is in fact better than its 
PIPC counterpart under the same average power threshold in 
terms of minimizing the PU capacity losses, which holds for 
the cases of both ergodic and outage capacities of the PU 
channel, with/without power control. To illustrate this result, 
we consider for simplicity the case without the PU link power 
control, i.e., the PU transmits with a constant power, Q, over 
all the fading states. Suppose that the PU link channel power 
gain is denoted by hp, and that from the SU transmitter to the 
PU receiver denoted by hgp. Next, consider the following two 
cases, where the interference power from the SU transmitter 
at the PU receiver, denoted by Igp = hgpPs, with ps denoting 
the SU transmit power, is fixed over all the fading states in 
Case I (corresponding to the case of PIPC), and is allowed to 
be variable in Case II (corresponding to the case of AIPC). 
For both cases, a constant interference power threshold F is set 



and is assumed to hold with equality, i.e., for Case I, Isp = F, 
for all the fading states, while for Case II, E(/ip^'') = F. 
Taking the PU ergodic capacity as an example, which can be 
expressed as (assuming unit-power receiver Gaussian noise): 

Let Cp^ and Cp^^ denote the values of Cp in Cases I and II, 
respectively. The following equalites/inequalities then hold 

where ( l22l i is due to the Jensen's inequality (see, e.g., IfTTl ) and 
the convexity of the function f{x) ~ log ^1 + where k 

is any positive constant and a; > 0. Thus, it follows that given 
the same average power of the interference, F, it is desirable 
for the PU to have the instantaneous interference power Igp 
fluctuate over fading states (Case II) rather than stay constant 
(Case I), to achieve a larger ergodic capacity. 

In general, the results in ifSTl reveal a new interference 
diversity phenomenon for SS-based CR networks, i.e., the 
randomized interference powers from the secondary network 
can be more advantageous over deterministic ones across 
different transmit dimensions over space, time, or frequency 
for minimizing the resulted primary network capacity losses. 
Further investigations are required on interference diversity 
driven DRA for CR or other spectrum sharing systems. 

C. Beyond Interference Temperature: Exploiting Primary Link 
Performance Margins 

So far, we have studied DRA for CR networks based on 
the IT constraints for protecting the PU transmissions. Given 
that the IT constraints in general conservatively lead to an 
upper bound on the PU capacity loss due to the interference 
from the SUs fT3l . 1521 . it would be possible to improve 
the spectrum sharing capacities for both the SUs and PUs 
over the IT-based methods if additional cognition on the PU 
transmissions is available at the CR transmitters. For example, 
by exploiting CSI of the PU links, the CRs could allocate 
transmit/interference powers more flexibly over the dimensions 
where the PU channels exhibit poor conditions, without de- 
grading too much the PU link performances. These PU "null" 
dimensions could come up in time, frequency, or space. Thus, 
the IT constraints could be replaced by the more relevant 
primary link performance margin constraints 1521 . I53II for 
the design of DRA in CR networks, in order to optimally 
exploit the available primary link performance margins to 
accommodate the interference from the SUs. Following this 
new paradigm, many new and challenging DRA problems can 
be formulated for CR networks. As an example, consider the 
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same setup with a pair of single-antenna PU and SU links 
over fading channels as in the previous subsection. Instead 
of applying the conventional AIPC: Fj{hspPs) < F, over the 
SU power allocation, we may apply the following PU ergodic 
capacity constraint ll52l 

Efl + -^^) >C', (23) 
V 1 + hspPs J 

where Cp is a given threshold for the minimum PU ergodic 
capacity. Note that the new constraint in (l23T l is more directly 
related to the PU transmission than the conventional AIPC. 
However, it can be verified that the constraint in (|23T l is non- 
convex over ps in general, thus resulting in more challenging 
SU power allocation problems than that with the convex AIPC. 
The optimal power allocation rules for the SU link subject to 
the AIPC vs. the newly introduced PU ergodic capacity con- 
straint given by (l23T l are compared in 1521 . where it is shown 
that the new constraint achieves notable rate improvements for 
both the PU and SU links over the conventional AIPC. 

V. Conclusions and Directions for Future Work 

Dynamic resource allocation (DRA) has become an es- 
sential building block in CR networks to exploit various 
cognitions over both the primary and secondary networks for 
CR transmit optimization subject to certain required primary 
protection. In this article, we have presented an extensive 
list of new, challenging, and unique problems for designing 
the optimal DRA in CR networks, and demonstrated the key 
role of various convex optimization techniques in solving the 
associated design problems. In addition to those open issues as 
highlighted in our previous discussions, other promising areas 
of practical and theoretical interests are discussed as follows, 
which open an avenue for future work. 

Robust Cognitive Beamforming: In our previous discussions 
on cognitive beamforming, we have observed that the knowl- 
edge of channels from each secondary transmitting terminal to 
all PUs is essential to the design optimization. However, since 
the primary and secondary networks usually belong to different 
operators, it is difficult for the PUs to feed back the required 
CSI to the CRs. As a result, the SU usually needs to rely on its 
own observations over the received signals from the primary 
terminals to extract the required CSI 1541 . Nevertheless, the 
estimated CSI on the SU-to-PU channels may contain errors, 
which should be taken into account for the design of practical 
CR systems. This motivates a new and challenging research 
direction on robust designs for cognitive beamforming to cope 
with imperfect CSI 1551 . l56l . More investigations on the 
robust cognitive beamforming designs for more general CR 
networks and CSI uncertainty models are appealing. 

Active Interference-Temperature Control: In this article, we 
have focused on the design of CR networks subject to the given 
interference-power constraints for protecting the PUs. We have 
also discussed some promising rules on how to optimally 
set the IT constraints in the CR network to achieve the best 
spectrum sharing throughput. These results lead to a new and 
universal design paradigm for interference management in CR 
or other related multiuser communication systems l57l . l58l . 



via appropriately setting the IT levels across the coexisting 
links. The active IT control approach to interference manage- 
ment for multiuser communication systems is relatively new, 
and more research endeavors are required along this direction. 
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